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COSIHC RAY SOURCES, ACCELERATION, 

AND PROPAGATION 

V.S.Ptuskin 

Institute of Terrestrial Magnetism, Ionosphere and Radio 
Wave Propagation, USSR Academy of Sciences, 142092 Moscow 


Review of selected papers on the theory of CR propagati- 
on and acceleration, presented in divisions OG;, 2, 5, 6, 7, 

8, and related problems, 

I. CR propagation in the Galaxy 

I. I, CR transfer in a turbulent medium. 

The high isotropy and a comparatively large age of galac- 
tic CR are explained by the effective interaction of relati- 
vistic particles with random and regular electromagnetic 
fields in interstellar medium. The kinetic theory of CR pro- 
pagation in the Galaxy is formulated similarly to the elabora- 
te theory of CR propagation in heliosphere (see the review by 
Quenby, 1984). The substantial difference between these theo- 
ries is explained by the necessity to take into account in 
some cases the collective effects due to a rather high density 
of relativistic particles. In particular, the kinetic CR stream 
instability and the hydrodynamic Parker instability is studied 
(see Cesarky, I960). 

The interaction of relativistic particles with an ensem- 
ble of given weak random magnetic fields is calculated by 
perturbation theory. The theory of CR transfer is considered 
to be basically completed for this case. (A new field of ac- 
tivities is suggested by Webb (8.3 - 8) in his paper on CR 
diffusion in relativistically moving plasma.) The main problem 
consists in poor information about the structure of the regu- 
lar and the random galactic magnetic fields. 

To calculate the diffusion coefficient of a particle 
with a gyroradius Th , it is necessary to know the spectrum 
of a random field SB* (k) in the resonant region of wave 
vectors k r<s ~ I/r M , The CR diffusion coefficient along a 
magnetic field is equal to 

T) - vr H _6, 


(I. I) 
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Remember that the gyroradius is 

V=~~ * 1.1 W 1 - (C/iC-teV) (6/3 - (v& c)S J - 2) 

The diffusion of CR with energies 10 - 10 eV is therefore 

determined by random fields with a size of inhomogeneties 
from 3 10 ,z - 3 10 20 cm (for E > 3 10'»E ®V, the gyroradius 

of particles exceeds the main turbulence scale in interstel- 
lar medium, ^ 100 pc, and the scattering becomes not very 
effective, D u E* ). 

The existence of the random magnetic field spectrum ne- 
cessary for CR diffusion has not been reliably established 
although it is not excluded by the available observations 
(Armstrong et al., 1981). At the same time, the simplest vers- 
ion of realization of a unique spectrum as an ensemble of li- 
near mhd waves seems rather doubtful because of the presence 
of a strong wave damping. 

In their papers Highdon (7.3 - 13) and Bykov and Toptygin 
(7.8 - 14) give a thorough theoretical analysis of turbulence 
formation processes in interstellar medium. The first-mention- 
ed paper deals with a small-scale turbulence due to evolution 
of isobaric entropy structures in interstellar plasma. Turbu- 
lence is shown to appear even on scales substantially smaller 
than the Coulomb mean free path. But the two-dimensional per- 
turbations (k6c« 0) arising in this case are inefficient for 
CR scattering. The second paper considers generation of secon- 
dary shocks appearing in the interaction of a primary shock, 
which is due to a supernova explosion, with interstellar 
clouds. The interstellar turbulence spectrum is determined. 

The corresponding value of the CR diffusion coefficient turn- 
ed out to be equal to D„ 2s 5 10 28 cm 2 /s. 

The problem of transverse diffusion of strongly magneti- 
zed particles in a stochastic large-scale magnetic field, which 
has not yet been strictly solved, is discussed in the paper by 
Ptuskin (7.2 - 16). 

Ginzburg and Ptuskin (7.2 - 15) present a microscopic 
calculation of the force acting on the CR gas on the side of 
background plasma and discuss the applicability limits of 
standard hydrodynamic equations, including the CR action. 

The relatively high CR pressure provides an important 
role played by relativistic particles in the formation of 
equilibrium distributions of gas and magnetic field in the 
galactic halo. This classical problem is considered anew by 
Dougherty et al. (7.2 - 18). The new element is an account 
of a finite pressure of mhd waves generated due to the stream 
instability of CR escaping from the Galaxy. The formation of 
an extended '’tail" in the spatial distribution of gas over 
the galactic plane is confirmed. 
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1.2. Semi empirical galactic models (E < 10 


15 


eV) 


Kinetic theory gives serious ground for using the diffu- 
sion approximation in tjie description of propagation of CR 
with energies up to 10* - 10'* eV in the Galaxy. Using the 
"microscopic" theory only, one cannot however strictly prove 
the diffusion character of motion and unambiguously determine 
the diffusion tensor and the velocity of OR convective trans- 
fer. Semi empirical models are of particular importance in this 
situation. They make it possible to classify and correlate nu- 
merous observational data, to explain the specificities of the 
composition, energy spectra, and anisotropy of different CR 
components, to find the CR composition in the sources. 

The diffusion galactic model is the most thoroughly deve- 
loped and on the whole explains well the relative observations. 
Its first basic version was proposed by Ginzburg and Syrovat- 
sky ( 1964 ) (for the modern version 3ee, for instance, the re- 
view by Ginzburg and Ptuskin (1985)* To simplify calculations, 
this model can in some cases (but not always!) be replaced by 
the leaky-box model. Many important problems remain, however, 
insufficiently investigated. 

At the Conference particular attention was given to the 
analysis of the role of particle acceleration in interstellar 
medium, to the study of stream instability of low-energy CR 
(Bretthorst and Margolis 7.2 - 9), to the clarification of 
consequences of strong interstellar gas density variations 
for CR transfer and fragmentation (Morfill et al. 7.2 - 4). 

The interest in acceleration in interstellar medium is 
connected with a great popularity of the scheme of diffusive 
shock acceleration by extended STJ remnants. The decisive argu- 
ment against the substantial CR acceleration in the course of 
their propagation and fragmentation in interstellar gas is the 
observed decrease with energy of the amount of secondary nuclei 
in the CR composition. 

The argument which was considered doubtless is objected 
in the papers by Lerche and Schlikeisei* (8.3 - 2, 7.2 - 8, 

8.3 - I;. According to formal calculations of these authors, 
in a continuous Fermi acceleration of CR in the entire Galaxy, 
the ratio of number densities of secondary to primary nuclei 
may decrease with energy. The crucial point of the model pro- 
posed by Lerche and Schlikeiser is the introduction of two 
different leakage times T^ - T^ p in the equations for concen- 
trations of secondary and primary nuclei on the basis that the 
spatial distributions of interstellar gas and CR sources do 
not coincide. This procedure does not seem to be correct. 

(Under some simplifying assumptions applied to this case the 
whole information on propagation of stable nuclei, both pri- 
mary and secondary, is contained in the particle distribution 
function with respect to pathlength x, G (x,r), here ? is 
the point of observation, see Ginzburg and Syrovatsky (1964). 
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Por the leaky-box model to be applicable, it is necessary and 
sufficient that the distribution function have the form 
G(x,j*)_^= exp (-x/x.d* )). In this case at each observation 
point r there exists only one value of x, (« ncT ), comfton for 
all nuclei, which should be used in transfer equations for CR 
concentration.) 

Using different methods, Giler et al. {8.3 - 4) and Cow- 
sik (8.3 - 7) have shown that a simultaneous CR acceleration 
and propagation in interstellar medium lead to sec/prim ratio 
increasing with energy (in some cases s/p -v const for E-*-*® ). 

Thus, acceleration in interstellar gas, which is accompa- 
nied by nuclear fragmentation, is excluded as the main process 
of CR acceleration. This does not mean, of course, that parti- 
cles accelerated in compact sources cannot undergo any additi- 
onal reacceleration. Such a scheme seems fairly probable (Sil- 
berberg et al. 1983, 8.3 - 5). 

Simon et al. {‘8.3 - 3) have studied a leaky-box model in 
which CR accelerated in their sources are then moderately re- 
accelerated during propagation. A model with reacceleration 
time approximately equal to leakage time from the Galaxy has 
been revealed to agree well with the data on the amount of se- 
condary boron nuclei in the interval of I - 100 GeV/n. In this 
case x t » 4g/cm 2 - (R C 6 GV) and X t co R“ 0 * 3 (R > 6 GV), 
whereas in the standard model without reacceleration X* = 

=» 8g/cm x and Xg <✓> R"* °‘ s (R > 6 GV). Thus, the real energy 
dependence of CR leakage time may be substantially weaker than 
in the standard leaky-box model (see Pig. I). This makes easier 
the interpretation of observations of CR anisotropy which point 
to a weak energy dependence of the leakage time for E < 10 /4 eV. 



Pig. I. Schematic drawing of the energy 

dependence of X * (E) as predicted 
by the standard leaky Box and the 
distributed reacceleration model. 
OG 8.3 - 3. 
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Note that reacceleration should not necessarily be asso- 
ciated with diffusive shock acceleration. If CR diffusion in 
interstellar medium is due to particle scattering by an iso- 
tropic mhd turbulence, the same scattering inevitably leads 
to statistical particle acceleration. It can be easily shown 
that the ratio of the particle leakage time from the Galaxy 
to the characteristic time of their acceleration is 


Tf/T, ^ - V~ (1.3) 

Here ^ is the Alfven velocity, J«fey(z) is the gas thick- 
ness across the galactic disk. With a standard choice of para- 


meters, we obtain T^ 
energy increases, 
Inessential. 


T a for GeV-energy particles. As the 
(E) decreases, and acceleration becomes 


1.3. Radioastronomical evidence 


Radi oastronomi cal observations make it possible, in 
principle, to establish the dimension of the CR confinement 
region in a Galaxy (the halo dimension) and also to find the 
diffusion coefficient and the speed of convective CR outflow 
(the galactic wind speed). It is difficult to interpret the 
radio-maps of our Galaxy since we are inside the radiating vo* 
lume. When observing other normal galaxies "from outside", it 
is often easier to distinguish the region of their halo. 
Examples are on-edge galaxies NGC 4631 and NGC 

891 in which clearly pronounced radio halos were revealed. 

Lerche and Schlikeiser (6.2 - I) and Cowsik and Sukumar 
(6.2 - 3) have constructed CR propagation models for NGC 46 31 
on the basis of radio continuum observations for frequencies 
from 327 MHz to 10.7 GHz. In the former paper it is stated 
that the observed dependences of the effective radio halo di- 
mension on the frequency and of the spectral radiation index 
on the height over the galactic plane are an unambiguous evi- 
dence of the existence of a large-scale galactic wind and a 
convective CR transfer in the galaxy NGC 4631. The latter pa- 
per shows, on the contrary, that a simple diffusion model 
without convection explains well the available radio data. 

1.4. Ultra high-energy CR 

As has already been mentioned, for particles with ener* 
gies E > 3 10 ,? a eV, scattering on inhomogeneities of the 
galactic magnetic field becomes inefficient and the diffusion 
coefficient increases rapidly with energy (D«4,r j| B£ /L* (STB)*- 
for r u > Lj. ). Diffusion gives turn to particle drift in 
an inhomogeneous regular magnetic field. At still higher 
energies, when the Lam or radius is comparable with the di- 
mension of the region occupied by the regular field, the 
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particle motion becomes simpler, and as the energy increases, 
it differs smaller and smaller from the free motion. For the 
field strength B - 2 10~ 6 G and for the dimension of the halo 
with an ordered field h = 5 Kpc, the condition r u < h holds 
for the energies E £ 10 ia ^ eVV 

An informative concise presentation of the origin of 
ultra high-energy CR is given by Hillas (1984). 

Most of the Authors believe at present that particles 
with energies 10 ,f - 10' a eV are of galactic origin, and for 
E>10* 3 eV there dominates the extragalactic component. The 
arguments in favour of this assertion are given in the paper 
by Efimov and Mikhailov (5.4 - 15) (this paper is a continuat- 
ion of the paper by Berezinsky and Mikhailov, 1983). The main 
attention is given to the weakest point of the galactic model 
- to the explanation of high CR isotropy. It is assumed that 
CR sources are distributed in the galactic disk, the regular 
azimuthal magnetic field in the disk B ^ = 2 - 3 10” 6 G, in 
the halo B y ^10“® G. The corresponding trajectory calculat- 
ions yield a good agreement with the observed values for the 
amplitude and phase of the I-st harmonic of CR anisotropy and 
for the southern excess of particles in the energy range from 
10 to 10 ‘9 eV. The apparent S/H asymmetry is probably connec 
ted in this case with the enhancement of particles from the 
general direction of the galactic plane (like in the model 
of galactic plane excess proposed by Wdowczyk and Wolfenda- 
le, 1984, 5.4 - II). . a 

For the energies E > 3 10 ,a eV the observed picture chan- 
ges drastically and there appears a particle flux from high 
northern galactic latitudes - roughly speaking, from the di- 
rection to the supercluster Virgo. This can be explained only 
by the action of extragalactic sources. 

The observed southern excess of particles with energies 
below 10 ,a eV can be explained, in principle, by a non-symmet- 
ric CR outflow across the galactic plane. In this case, the 
particle concentration in the region P in Fig. 2 must be lower 
than in the external region Q. The Hillas interpretation 



Fig. 2. Trajectories of positively charged particles in 
azimuthal galactic magnetic field. OG 5.4 - 9. 

(1984) suggests the presence of a local CR gradient in this 
energy range (~ 15 % Kpc ) increasing toward the Orion regi- 
on. Sommers and Elbert (5.4 - 9) believe that in the region P, 
at a distance of several Kpc from the solar system there oc— 
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curs a reverse of the azimuthal field which is ^ 

a rapid evacuation of CR and a local decrease of their inten- 

Sity *Karakula and Tkaczyk (5.4 - 8) have calculated the tra- 
jectories of ultra high-energy particles and determined tran- 
sparency of the Galaxy for extragalactic protons. The regular 
galactic magnetic field is assumed to be contained in the disk 
laU 0.4 Kpe. In constructing a concrete model of CR origin, 
the authors assume the isotropic flux from galactic sources 
to be summed up with the directed flux from the source >*0 \g 
Virgo Cluster (these fluxes are equal to each other for io 
eV). 


2. Antiprotons 

An understanding of the large flux of antiprotons found 
in CR remains an intriguing problem. The available observat- 
ions (unfortunately, they have been little enriched during 
this Conference) can be triefly resumed as follows. The ob- 
served integrated antiproton flux was found to exceed signi- 
ficantly the flux of secondary antiprotons calculated from 
the standard model of CR propagation in interstellar medium. 
(See, for example, Stephens 1981). For energies 0.1 - 10 GeV, 
the observed total flux of antiprotons 4—10 times exceeds the 
flux expected in the leaky-box model for - 5 g/ cm . -he 
observed energy spectrum of anti protons is also unexpectedly 
different from the production spectrum of secondary antipro- 
tons generated in p— p collisions. It is greatly enriched with 
low-energy particles and on the whole is similar to the CR 
proton component spectrum (in contrast with the production 
spectrum of secondary antiprotons which falls sharply for the 
energies E«£2 GeV). 

The explanations proposed may conditionally be divided 
into ’’exotic” and ”non-exotic” . , _ 

Here is the list of ’’exotic” explanations exploiting the 
new physical principles: quantum evaporation of mini black 
holes Kiraly et al. (1984), n-n oscillations Sawada et al. 
(1981), Sivaram and Krishnan (1982), primary extragalactic 
origin Stecker et al. (1981), Stecker and Wolfendale (6.1 - 
8), photino annihilation in the galactic halo Silk and Sred- 
nicki (1984), Stecker et al. (6.1 - 9). 

The photino hypothesis is based on the assumption that 
in the Nature there exist stable massive photinos. These 
particles originate in the early Big-Bang. Photino is assumed 
to make up the missing mass in the galactic halo and to pro- 
bide the matter density in the Universe which is close to the 
critical one SI * I. In the paper submitted to the Conference 
Stecker et al. investigate photino annihilation in the galac- 
tic halo and calculate the spectrum of antiprotona appearing 
in this process. An intriguing fit is obtained to all the 
existing data on antiprotons for a photino mass mcrss15 GeV. 

The cut-off of the spectrum for E * Vwy is predicted. (One 
should remember that it is still unknown in what concrete 
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form, if at all, the very supersymmetry principle is realized 
in the real world) . 

Non-exotic hypotheses do not require new physical mecha- 
nisms for explanation of a high p flux. However, they require 
reconsideration or specification of the old habitual galactic 
models of CR propagation (see the Discussion in the papers by 
Ginzburg and Ptuskin, 1984; Lagage and Cesarsky, 1985). 

It is possible, in particular, that part of CR sources 
(or all the sources -at a certain stage of their evolution) 
are surrounded by a thick layer of matter in which secondary 
p are produced. Heavy nuclei either do not escape from such 
objects due to a strong fragmentation or are not accelerated 
in them at all. One of the possible realizations of such a 
model has been considered by Mauger and Stephens (1983). 

Primary CR are assumed to accelerate in SN explosions in den- 
se clouds and to be confined there for several thousand years, 
traversing the thickness of about 50 g/cm^. If -v 30 % of nucle- 
ons observed in CR come from such sources, one can satisfacto- 
rily explain the observations of p. Now Stephens shows (6.1 - 
7, 6.2 - 9, 2.5 - 3) that the fluxes of secondary positrons 
and gamma-quanta expected in this model do not contradict the 
observational data. More favourable here is the version in 
which the time of CR leakage from the Galaxy is almost inde- 
pendent of the energy. If the dependence is strong (T^ooE " 0,6 ) , 
the flat CR spectrum in the sources leads to a too rigid spec- 
trum of gamma-rays and to a too large flux of secondary posi- 
trons. 

Derraer and Ramaty (6.1 - 4) have presented a model in 
which low-energy antiprotons appear as secondary in p - p in- 
teraction in a relativistic plasma with a temperature kT ~ 

0,2 m p c z . It is assumed that the appropriate conditions may 
exist in the vicinity of a neutron star or a black hole. 

Various versions of secondary anti pro ton generation in 
dense gas clouds in models with a nonuniform CR propagation 
are developed by Tan (5.4 - 13, 6.1 - 6, 6.2 - 7), Dogiel et 
al. (8.2 - 17), Morfill et al. (7.2 - 4). 

The model which would explain the high antiproton flux 
in the CR composition has not yet been finally chosen. 

The correctness of calculations of the expected fluxes 
of secondary p has been verified by Bowen and Moats (6.1 - 3). 
The parameters «( = 0.333 - the probability of np charge ex- 
change, and £/ 2» 0.45 - the average elasticity, have been 
determined by measuring the proton spectrum at mountain alti- 
tude. The calculated spectrum of secondary p in the atmosphe- 
re satisfactorily agrees with the measurements by Bowen et al. 
(1983). The contribution of atmospheric p in balloon measure- 
ments by Golden et al. (1979) and Buffington et al. (1981) 
has been shown not to exceed about 10 %, 

3. Shock acceleration: theory and application to the 
CR origin problem 

The diffusive shock acceleration mechanism remains the 
moat popular with theoreticians engaged In the problem of 
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CR acceleration (see the reviews by Axford, 1981, Drury, 1983, 
Volk, 1984). 

3.1. Test particle approximation 

An acceleration of a fast test particle diffusing near a 
shock front is a version of a first-order Fermi acceleration. 
It is due to a repeated crossing of the shock front in a ran- 
dom particle walk and to an energy gain in front collisions 
with scattering centres embedded to the background plasma. A 
formal solution of the problem can be obtained using the equa- 
tions for the test particle distribution function which desc- 
ribes spatial diffusion, convective transfer, and a regular 
alteration of particle energy in a nonuniform flux: 

|| + J - vt? -g |i = o, (3.1.) 

The particle density is equal to N(p) dp = 4irp*f(p) dp. 

For a plane front the velocity profile u(x) is approxi- 
mated by a step function (see Fig. 3). If particle distributi- 
on in a nonperturbed medium in an upstream region has the form 
f .go S' (p - p 0 )« the stationary spectrum of accelerated par- 
ticles in a downstream region has a power-law form: 

3r 

©(p-poV (3.2) 

where r * u^/u^ is a compression ratio in a shock. 



Fig. 3. Spatial distribution of accelerated teat 
particles in the shock frame 


For an extremely strong shock wave without radiation, 
which propagates in a gas with the adiabatic index jo.* 5/3* 
the quantity r » 4 and, therefore, the spectrum of accelerat- 
ed particles is N(p) co p"*. This is close to the expected CR 
spectrum in the sources. 

The characteristic time of particle acceleration (i.e. 
the time of its energy variation by a factor of e ) under 
shock acceleration is equal to 



3- (D t. 

* 



(3.3) 


Note that the estimate of the characteristic accelerat- 
ion time t.^D/uHs valid both for the diffusive shock acce- 
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leration and for the statistical Fermi acceleration. In the 
latter case u haa the meaning of the turbulent motion velo- 
city. 

If a particle bears additional energy losses in a medium 
(with a characteristic time t £ ), then for acceleration it ia 
necessary that the condition t a c< should be fulfilled. In 
particular, for a spherical shock wave with a radius R, the 
role of losses, is played by a diffusion particle escape for 
the time t < i~R < Vl), and therefore particles are accelerated 
only under the condition (the numerical factor is omitted 
here and it is assumed that B^>, B^) 

RixJ > i. (3.4) 


These assumptions are confirmed, in particular, by the nume- 
rical calculations made by Ko and Jokipii (8.2 - 2). 

By the terminology accepted in the theory of particle 
propagation in a solar wind, the inequality (3.4) implies a 
strong GR modulation. 

Since the diffusion coefficient usually increases with 
energy, the condition (3.4) imposes limitation on the maximum 
possible accelerated particle energy E m . Under scattering 
in a magnetic field D w > r H v , accordingly. 


j£ t\r 

max 


%eB 


(D t ~r H ir, u~c). 


(3.5) 


An impressive example of a possible realization of the 
mechanism of CR shock acceleration in the Galaxy is given by 
Jokipii and Morfill (8.1 -8). They suggest that CR particles 
with energies E ^ 10' 5 eV are accelerated at a termination 
shock of galactic wind. If the wind in the Galaxy does actual- 
ly exist (which can be strongly doubted, see, for instance, 
Habe , Ikeuchi, 1980) and if the large-scale magnetic galactic 
field has a structure similar to the Parker spiral in inter- 
planetary space, then B ~ B c (R©/R) for R > R 0 , where B c = 

3 10~®G is a field at a distance R©= 10 Kpc. It is assumed 
that the distance to the shock is R ~ 100 - 200 Kpc, the wind 
velocity u ^ = 300 - 500 - km/ 3. 

Formula (3.5) does not work in this case since in an or- 
dered, strongly twisted spiral magnetic field, the diffusion 
coefficient in the radial direction is smaller 


X>ejj = <L Z I> U + (3.6) 

where = Rol^vC RHjal) (“Ujai = 250 km/ s is the speed of 

Galaxy rotation), D A ~ r£cVD„ is the diffusion coefficient 
across a regular magnetic field. 

The minimum value of D e t$ is reached for r H c/oL-v 
20 r^c. In this case, instead of (3.5) we have 
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The maximum energy density of accelerated particles can 
be estimated by assuming that the whole power, of the galac- 
tic wind (Q~5 10*° erg/s, according to Jokipii and Morfill) 
goes to accelerated particles 

*,/«* <3.e) 

Here the time of CR escape from the galactic wind region is 

T^R/liv ~ 3. 10** (R/ 100%} years. 

The CR intensity modulation effects lead, in fact, to a 
substantial lowering of the uftr value. In this case, par- 
ticles with E £ 10'® eV do not reach the observer at all. We 
should recall that the values observed near the Earth are 
-W cr (>10 3 eV) 5c 10 - "' 2 erg/cm a and %«/<;,.(> 10'® eV) 

5 10“' 6 erg/ cm 3 . 

Thus, CR acceleration on the galactic wind boundaries 
could, in principle, be a noticeable source of particles 
with energies E > 10*® eV. It remains unclear how the main 
energy output can be provided just in the region E > 10 ,y eV 
(with the particle spectrum N(E)dE eoE‘ M dE) and how for 
E < 10'® eV a smooth conjugation with the CR spectrum (JT(E)co 
E“ 2 * 7 ) generated, according to Jokipii and Morfill, in galac- 
tic 3N remnants can be obtained. Rote that reacceleration of 
CR supplied by SR is inessential in this case because due to 
the geometrical factor it gives not more than (R® /R)*^ 
of the CR concentration in the galactic disk. It is necessary 
to consider particle acceleration directly from the thermal 
plasma on the galactic wind boundaries. 

SR remnants are regarded as "classical" astronomical 
objects, in which the action of diffusive shock acceleration 
is possible. A direct interpretation of radio data from SN 
remnants with an age exceeding 100 years is ambiguous. The 
hypothesis concerning diffusive shock acceleration is consi- 
stent with observations, but is not rigorously proved (Beck 
et al. 8.1 - 10, Bogdan et al . 8.1 - II, Lawson et al. 6.2 - 
4 ). 

Old remnants of the SR Loop I and III are itudied in the 
paper by Lawson et al. Diffusive shock acceleration of elect*-' 
rons up to energies E<^10 GeV has been revealed. The CR diffu- 
sion coefficient in the galactic disk has been determined to 
be D ~ 10*^ cm 2 /s. 

Of importance is finding the fraction of SN explosion 
energy which can be transferred to CR. This problem has been 
solved numerically in the test particle approximation in the 
paper Jjy Dorfi and Drury (8.1 - 9). With the explosion energy 
of 10 5 ^erg and with the diffusion coefficient equal to 10 AT 
cm^/s, the effective acceleration of background relativistic 
particles starts at the moment t * 6 10 sec after the 
explosion, and by the moment t * 2 10** sec about 12 % of the 
initial SN energy is transferred to high-energy particles. 
Such an efficiency is, in principle, sufficient to replenish 
the observed CR energy density. 
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3.2. Nonlinear CR shocks 


The high efficiency of shock acceleration and a compara- 
tively large energy density of relativistic particles in cos- 
mic plasma lead to the necessity to study the back reaction 
of accelerated particles to thermal plasma. (In this case the 
quantities u and D in equation (3.1) cannot be regarded as ex- 
ternally given parameters). This nonlinear problem turns out 
to be very complicated. 

If CR are assumed to be a separate relativistic gas com- 
ponent with a pressure P cr , an internal energy density v cr 
and an effective average diffusion coefficient D 


Pcr= > 


w 


cr 


oo 

= 4t r u 


Z r 


Up), 


(3.9) 


c>0 cx? 

sun 

then in the double-fluid hydrodynamic approximation the equa- 
tions of one -dimensional motion of a medium, with an account 
of CR action, have the form 


2JL 9 


(?“)= 0, 



v'is 

*>X 



(3.10) 


^ VX U > + 1 ^7 ° cr -rf x p cr . 

In the study of self-consistent shock structure, one 
seeks for steady solutions which would describe the transition 
from the state given for x =-oo (upstream) to the state uni- 
form for x = +oo (downstream), see the reviews cited at the 
beginning of Sec. 3. 

A more complicated problem, which takes into account the 
finite energy of CR-scattering mhd waves, has been considered 
in the papers by Volk and Me Kenzie (1982), Volk et al. (1904). 
These waves are assumed to be generated due to stream insta- 
bility of relativistic particles before the shock front. The 
set of equations (3.10) must then be modified in the respecti- 
ve way. 


An account of CR pressure changes the profile of a hydro- 
dynamic flow near the shock front, see Pig. 4. There appears a 
region of smooth variation of upstream velocity (precursor), 
and in the general case there remains a step-like jump of ve- 
locity (subshock) determined by gas viscosity. This situation 
is similar to isothermal jumps for strong shocks in media 
with a thermal conductivity coefficient substantially exceed- 
ing the viscosity coefficient in the usual hydrodynamics 
(Landau and Lifshitz, 1959). 
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Pig. 4. A stationary modified shock 


Unfortunately, one had not yet succeeded in obtaining a 
single-valued steady-state solution to the problem (more pre- 
cisely, in obtaining generalized Rankine-Hugoniot conditions 
for arbitrary values of N = P cr (-00 )/(P t r (-<*>) + PjC- 00 )) 
and the Mach number of the incoming flow M£ = 9 u£ * \ + 

+ V P )‘% It is shown, in particular, that for = 5/3, 

= r 4/3 in the general case there may exist three soluti- 
ons 'in the downstream region for certain upstream conditions. 
Such an unambiguity is caused by divergency of the quantity 
P er with the momentum p for flat spectra f(p) to p-^ ( 0), 

See (3.9). Then the definition (3*9) also loses sense for 
D (if D(p) increases with p). Indeed, let the diffusion coef- 
ficient increase with energy and the transition region from, 
u ( — 00 ) to u(+ 00 ) have the finite dimension AX . For partic- 
les with momenta p » p* , where the value p* is determin- 
ed from the condition u t .A. x/D(p*) = I, acceleration proceeds 
in the same wav as on a step-like profile of the velocity 
u(x) * u t 0(-x/ + Ug0(x), and the spectrum of high-energy 
particles has the form (3.2). The degree of compression for 
a strong shock exceeds here r = 4 since the presence of rela- 
tivistic particles in the downstream region softens the equat- 
ion of state ( Ycr < XV ) • Ia ' thia case “ he spectrum f(p)co p“'’ + £ 
and P fr = 00 . This implies that a steady-state aolution with 

a limited total width of the transition region is impossible. 
Additional difficulties are caused by the functional dependen- 
ce y er (f(p)). 

It has become clear that in general either the maximum par- 
ticle energy must be restricted by introducing additional loss 
proccesses into the problem or the acceleration must be treated 
as a time-dependent problem (Drury, 1984), 

The time-dependent structure of a nonlinear CR shock is 
numerically investigated in the paper by Dorfi (8.1 - 3). For 
not too large a value M* =4.4 and for N = 0,5, the transit- 
ion from a gas dominated shock to a shock modified by CR ta- 
kes up the time of ( 30 - 40) D/u z . For a high Mach number 
17 = 10 and for a small value N = 0 . 05 , there appears no back 
effect of the CR pressure on the motion of background plasma 
up to a time t ~ 10* D/u^r 

Beck and Drury (8.1 - 4) analytically investigate the. 
time -dependent problem. They construct a selfsimilar solution 
which depends on the similarity variable 7 ■ x/t. Selfsimilar 
shock structures are shown to contain always a subshock. The 
solution to the problem for shocks with a large Mach number 
has not been obtained. 
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The problem of seeking for an unambiguous solution, 
which has not yet been finally solved, is in a sense only of 
academic interest . Volk et al. (1984; have shown that in the 
framework of three-fluid hydrodynamics in the strong wave 
damping approximation the solution is unique in a wide range 
of reasonable astrophysical parameters. 

The study of the dynamics of SU remnants with an account 
of the action of accelerated CR, which is of great importance 
for astrophysics, has been started in the paper by Volk et al. 
(8.1 - 12), First results have been obtained on the time evo- 
lution of the kinetic remnant energy, the thermal energy of 
a heated gas, and the total CR energy. It has been shown that 
10 - 50# of the initial energy of an explosion can be trans- 
ferred into the CR energy. 

Another, not hydrodynamic approach to nonlinear CR shocks 
is developed in the papers by Eichler and Ellison (8.1 - 6) 
and Be rejthko et al. (8.1 - 13) (see also Eichler (1979, 1984), 
Krymsky (1981), (1983), Ellison and Eichler (1984)). This 
approach consists in a unified description of thermal and ac- 
celerated particles, i.e. along with acceleration one consi- 
ders thermal particle injection. Particle acceleration is an 
inherent part of the very process of Ihock formation in plas- 
ma. It is assumed that there exists a maximum particle energy 
E m<U( above which they cannot be confined near the shock. 
Therefore, there exists a continuous energy flux from a sy- 
stem x 

= l N W (3.9) 

— oO 

(j = 2 for nonrelativistic particles and j = I for ultra re- 
lativistic particles). The diffusion coefficient D(p) is as- 
sumed to increase with momentum. 

For fast particles the distribution function obeys equa- 
tion (3.1). Using the mass and momentum flux conservation con- 
ditions for given upstream values, one finds a stationary 
solution of equation (3.1). The final solution is obtained by 
matching the low-energy spectrum region to the thermal back- 
ground distribution. 

An alternative procedure consists in the numerical solu- 
tion of the kinetic equation with a simplified collision in- 
tegral in the entire particle momentum region and in the defi- 
nition in this way of the complete distribution function. 

In the paper by Eichler arid Ellison (8.1 - 6) in the des- 
cribed scheme, the final velocity of waves generated by an ac- 
celerated particle stream i3 taken into account for the first 
time. Figures 5,6 present the obtained spectrum of accelerat- 
ed particles for the case of acceleration. .by SIT remnants in 
interstellar medium. The values E^^ =» 10 ,5 eV, B * 3 10"®^ 5, 
n = lO-^crn" 3 , u± = 10 3 km/s were used. 

It is remarkable that for a wide range of reasonable as- 
trophysical parameters the spectrum turned out to be close to 
E(E k ) o) . Typical is also the presence of two maxima in 
the function EA N(E^ ) - for thermal energies and for E . 
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For the spectra shown in Fig 5 the kinetic energy flux going 
into relativistic particles makes up 72 - 26 %, 

Such a calculation is in a good agreement with observa- 
tions for shocks in interplanetary medium (Ellison and Eich- 
ler, ’M984), but the problem as a whole cannot yet evidently 
be considered as finally solved. Suffice it to say that the 
indicated scheme must include, in particular, the theory of 
collisionless shocks in a usual plasma. This theory has been 
intensively developed for already about 30 years. 



LOG ENERGY (eV) 

Fig. 5 Postshock partial pressure va. 
energy. An energy cutoff = 10 ^eV has been 
used. Each label represents a family of 
curves with a particular u i/u(GeV) .06 



MACH NUMBER, M-, 

Fig. 6 Labels same as Fig. 5 
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It ia of importance to note, however, that the difficul- 
ties in the theory of nonlinear CR shocks are due to the high 
efficiency of the diffusive shock acceleration and, in fact^ 
we try to understand what it is limited to. In any case, the 
efficiency of the conversion of the energy of an ordered hyd- 
rodynamic stream into the energy of relativistic particle's 
typical for the galactic model of CR origin ( /v 10 %) does 
not seem to be excessive. 

There are no essential difficulties with explanation of 
the observed power-law energy CR spectrum. The problem con- 
sists in obtaining a sufficiently large 10* 7 - 10* 9 for 

galactic CR souces. So, for Instance, the value io ,s eV 

presented in Fig. 5 has been not calculated, but postulated. 
For the values B = 3 10 G, n= 10~ 2 cm -5 used by the auth- 
ors, an SN explosion with a total energy W = 10 51 erg may. 
in accordance with (3.5), actually give E ^ =: 3 10™ eV. In 
the standard model of CR propagation in interstellar medium, 
the v^lue of D„/vr H v entering ( 3 . 5 ) seems to be strongly 
understimated • The value of E should rather be lowered by 
several orders of magnitude (for more details see Volk, 1981). 

The effect somewhat increasing E is generation of a 
strong random magnetic field near the remnant due to the stre- 
am instability of accelerated particles. This process cannot 
be investigated within the theory of weak turbulence, and its 
analysis remains a challenge for theoreticians (see Drury 
1933, Volk 1984). The possibility of hightening E ^ (may be 
up to 10 eV) under CR acceleration in the reverse ^shock 
inside an SNR is pointed out by Volk et al. ( 8.1 - 12 ). 

3.3. New problems 

In their paper ( 8,1 - 5 ) Dorfi and Drury have revealed a 
new interesting effect which accompanies CR diffusive shock 
acceleration. It has been shown that in the region of precur- 
sor there develops an instability of compressions! disturban- 
ces of the medium, due to the CR gradient, for the wavelengths 
* L “ 30/ v is the mean free path of a particle, 

h 7 D ' u „ is the dimension of the precursor). The instability 
arises for L ■< | I + S' O-n D)/ F(ln )|:'* D/v^. and develops 
during the time “t** L/ v S o M ^ . For I, the time 'f is 

small as compared with the time of convectional outflow of 

Q in ’r n q n ci -P v» rvrvi -r / tr . . . . 


' sr ~ - **--»** V/UiivcubiUimi UUfclXOW 01 

disturbances from the prcursor region L/v s „ M, . Slight densi- 
ty di sturbances in the upstream med5_um must be substantially 
amplified in the precursor region# 

Zank and Me henzie (8*1 - ?) confirm the existence of 
instability. They also investigate instability ©f a 
CR shock with respect to long-wave disturbances ( X>> L). In 
the presence of a relativistic CR gas, strong shocks always 
prove to be unstable. J 

The effect of these instabilities on the shock structure 
has not yet been investigated. One may expect a strong sto- 
chastization of the motion of the medium before a shock, the 
appearance of secondary shocks, an additional widening of the 



23 1 


front and even destruction of its plane structure. 

Webb ( 0,1 - I) was the first to study the structure of 
relativistic CR shocks. He points Outlhai, the effective CR 
acceleration proceeds for Ui 4 i c Joc r “ ^ ~ 8.58c. In nigh 

velocity streams (u 1 > c J j cr - 1 ) CR cool down. 

4 . High-energy particles in various astronomical 

objects 

4.1. Active galactic nuclei 


A.GN and quasars are possibly the most powerful CR sour- 
ces. It is believed now that the particles of ultra high 
energies E > 10 19 eV observed near the Earth have been acce- 
lerated just in AG N. To the presence of relativistic partic- 
les there testifies nonthermal AGN radiation (see Rees 1984), 
but concrete mechanisms of particle acceleration remain un— 


d P QX* • 

Kazanas and Ellison (8.1 - 7) consider diffusive shock 
acceleration of protons in accreting matter near a black hole 
of mass M 3 = M/(10 3 M 0 ) . Accelerated relativistic protons 

provide the pressure to support the standing shock. They un- 
dergo inelastic nuclear collisions and generate secondary 
* , V , which are responsible for the observed^ radia- 


^ f Q t V f JU CAJ ^ J. V. X v 

tion. For a nearly 100 % conversion of the hydrodynamic mo- 



is the Eddington luminocity, x t » R/R* is the radius of 
spherical shock in 3 chwarzschild radii. The model provides a 
natural explanation of the observed L <-o M correlation for qua- 
sars and galactic nuclei. In this case 10 ^ Xjl^. 200. For the 
case NGC 4151 (M = 3 10 7 M 0 , L a 10 ™ erg/s) L « 2.5 10 L E , 


x ^ = 140 . 


4 . 2 . Neutron stars in close binary systems 

ir 

Ultra high-energy (> 10 eV) - rays (or/ and unidentifi- 
ed neutral particles) have been observed from Cyg X - 3 and 
possibly from the X-ray binary sources LMC X-4, Vela X-I. 
Taking the measurements at face value, we have indication of 
a striking efficiency of CR acceleration in the compact bina- 
ry systems containing n-stars. It is sufficient to have CygX— 

3 alone to maintain the present flux' of the galactic CR above 
10*6 eV. The extragalacti c source LMC X-4 (Ey > 10 ,b eV) has 
(?) a luminosity more than 20 times the one of Cyg X-3. 

Taking the Haverah Park flux for Cyg X-3 Fy ( > 2 10 eV) 
= I.I 10 “' r erg/cm*s, Hillas (5.4 - 7) obtains the following 
estimate for proton luminosity of the source 


Lll0*-I0'\v)=£ * 6 * 3 » 50 >U<rr*xl.Ho' l< ’*& IMO^eq/* 1 ^* I > 
4ir 

Here we have the factors taking into account the appearance 
of particles in a solid angle SI * *he efficiency of energy 
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conversion p •*+ y , absorption in rout due to e“ e* pair crea- 
tion on relic radiation, pulse du±y ratio of the 4.8 h period 
of the system. The distance to Cyg X-3 is taken to be «qual 
to r = 12 Kpe; 

The proton luminosity being so high, Cyg X-3 must also 
be a source of secondary neutrinos (see Beresinsky 2.1-7 
Brecher and Chanmugam 2.2 - 5 , Gaisser and Stanev 1985). 

_ The hypothesis on propagation of neutral radiation from 
Cyg X-3 in the form of photino (V. J.S tenger) requires an im- 
probably high proton luminosity L p = 2 10^ («Q,/10 -2 -) era/s 
(Beresinsky 2.1 - 7). 

As the energy source of accelerated particles in a bi- 
nary system consisting of a normal star and a magnetized ne- 
utron star, one considers two main possibilities: I) accre- 
tion from a normal star onto an n-star; 2) pulsar action of 
the n-star itself. 

In any case, proton acceleration proceeds in the vicini- 
ty of an. n-star. The primary proton spectrum has a power-law 
form or is a monoenergetic beam with an energy E ~ 10*' eV. 

In the latter case, the power-law spectrum of radiation and of 
secondary particles is caused by a cascade in the target re- 
gion (Hillas 1984, 5.4 - 7). 

Secondary UEE % -rays appear as a result of pp collisions 
in an accreting gas outside the acceleration re gL on or in the 
atmosphere of the companion star. (Aharanian et al. (2.6 - 13 ) 
notice that at a rather high density of background low-fre- 
quency radiation in the source, a substantial contribution 
Is made also by photomeson X -ray generation). 

In the source Cyg X-3 in different periods of observa- 
tions in the UHE range, % -pulses had two different values 
for the phase of the 4 . 8 -hour binary orbit: = 0.25 and 

^ According to Hillas (5.4 - 7), in the periods of 

'f'vO.b the target for the production of secondary X -photons 
is a wake which occurs if accretion appears from a stellar 

w 9ee 3 i m il ar scheme is discussed by Protheroe 
and Clay (1985) for LIvIC X-4) . 

The dynamics of the interaction between an accelerated 
particle beam and the atmosphere of a companion star is con- 
sidered by Beresinsky (2.1 - 7). He proposes a heating model 
in whi ch radiation in the phase 'Pa. 0,2 is not accompanied by 
the appearance of a symmetric pulse 'P = 0.8. 


v\ n-2 

. -'P 1 Figure 7. Supposed geometry 
. for Cygnus X-3. 06 SA -7. 


y ** protons 


<>■«/ 
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To choose a concrete model of OR acceleration in binary 
systems, it is important to establish whether or not Gyg X-3 
is a unique UHE X -ray source. 

So, for a pulsar model (Eichler and Vestrand 1984, Be re- 
al nsky 1979, 9.1 - 7) a very young pulsar with a rotation 
speed of 1X3*- 10 3 times per sec. is needed. Such a version 
is excluded for the sources LMC X-4 (P = 13*5 a), Her X-I 
(P si 1.24 a), 4U OII5 + 63 (P = 3.61 s), but is not excluded 
for Cyg X-3 (the value of P in unknown). 

In accretion models of CR acceleration two schemes are 
considered - the dynamo model in an accretion disk and shock 
acceleration in an accretion flow. 

Brecher and Chanmugam (1985; 7.2 - 5) proposed a unipo- 
lar induction model for CR acceleration in accretion binary 
systems, such as Cyg X-3, Her X-I, Vela X-I, LMC X-4. 

Details of the electrodynamics of such systems are not 
yet clear. Particle acceleration is assumed to be caused in 
the end by a very high potential drop across the accretion 
disk between the Alfven radius and the external .edge of the 
disk. The scaling law L (L = GMM/R is the 

total accretion luminosity of an n-star of mass M and radius 
R, B is a magnetic field on its surface) is obtained for the 
maximum possible energy of an accelerated particle. The 
maximum luminosity in relativistic particles is estimated to 
be Lpy^^PcEtywnx /e* (e is the particle charge). Thus, for a 
given magnetic field strength the higher the accretion rate, 
the larger the maximum CR particle energy and the higher the 
total particle luminosity. For B = 10' 3 G the energy E v>)ax ~ _ 5 
10*' eV is reached for a very high accretion rate $=10~ 6 - KT’ 
M® yr 1-1 , iie. L<^10 w - 10 . erg/s. In this case to relati- 

vistic particles there goes L p wx^^l 0 33 erg/s. Most of the 
energy must be released in the form of a jet. In the unipolar 
induction model comparatively weak fields B/v 5 10® G are 
preferable. In this case E mx = 10' 7 eV, L p mftx = 10 39 erg/s ~L. 
The value B = 5 10® G is not excluded for Cyg X-3. (the field 
is unknown), but is excluded for Her X-I (3 = 4 10 lz G) and 
4U0015 + 63 (3 = 10 a G). 

The model of shock acceleration of particles in applica- 
tion to the source Her X-I is considered by Eichler and Vest- 
rand (2.2 - 8). For an explanation of % -radiation outbursts 
observed for energies F > 10' 2, eV and E > 10 ^ eV, they sug- 
gest a diffusive shock acceleration in an accretion column at 
a distance R ~ 30 radii of the n-star (at smaller distances 
acceleration is hampered by synchrotron losses). In this case 
E ~ 1 0 P me* (R/10 6 cm)- 1 ' 6 . The observed y -rays are 
formed in the interaction between protons and the surrounding 
accretion disk. The model predicts } ( -ray outbursts at the 
onset and decline of the high-intensity X-ray state. At this 
time our line of sight is grazing the accretion disk. Radia- 
tion of -rays with an energy E > eV is impossible in 

this model (see, however, Kazanas and Ellison, 1985). 

The existence of such UHE )f -ray sources as Cyg X-3, 
possibly leads to the formation of a galactic X-ray halo 
CRana et al . 2.2 - 6). Under yy collisions with relic photons, 
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UHF) X -rays produce e~ e + pairs, and these will in turn pe- 
magneticfield^ 3 * * ^ 3ynchrotron radiation in the galactic 

4.3. Stars 

Normal stars may he injectors of fast particles which 
are further accelerated up to high energies. Bogdan and 

A Sia r I 8 ** “ 9) consider resonant statistical acce- 

4^ in f 1 aring atars* The maximum particle 

momentum is limited m this case to the energy losses due to 
synchrotron radiation and in this model makes up p ^ 1,5 10® 
e//c for the magnetic field of a star B - 100 G. 

4.4* SN remnants 

Most of the papers on CR acceleration in SNR have alrea- 
dy been discussed in Sec. 3. 

In some SNR (Crab), relativistic particles are supplied 
by a central source - a pulsar. The remnant G 59.7 - 0.3 (the 
distance - 19 Kpc, the dimension - 1.8 pc) may also appear to 

8 vTi b + f thls type * resolution maps obtained with 

the VLA show two spectrally distinct components - a flat-spec- 
trum core surrounded by a shell. Koch-Miramond et al. (6.2 - 
1 *?Y Q presented the data of Exosat observations of X-rav 
radiation of the remnant in the 2-10 KeV range. The search 
for pulsations with.a period between 32 ms and 10 a in the 
radiation of the object has not yielded a positive result* 
Assuming radio and X-ray radiation to be, synchrotron, they ob- 
tained the values of the field B = 2 lO^G and of the total 
electron energy W -1.6 10«erg. A continuous energy inject- 
ion from the central source follows from the small synchro- 
tron electron lifetime as compared with the SNR age. 

4.5. Molecular clouds. 

the paper by Dogiel et al. (8.2 - 17) giant molecular 
clouds are interpreted not as passive targets for relativis- 
tic particles, but as active objects which accel&rate CR up 
to /w io GeV. The energy sources are hydrodynamic motions of 
a neutral gas which generate a chaotic magnetic field. The 
galactic CR that penetrate from without undergo a statistic- 
al acceleration, as a result of which the CR energy density 
increases by about an order of magnitude. This explains a 
hightened )( -luminosity of some molecular clouds (the expect- 
5 -ray spectrum (E) c* E* 1 * for E > 100 MeV). An inten- 
allied production of secondary particles in clouds could 
explain also an anomalous amount of p, d, 3 He, e* in CR. 

1 ra03 ' t sincerely thank the organisers for their risk to 

invite me as a rapporteur at 19 ICRC. Very useful discussions 

with many colleagues are gratefully acknowledged. I am thank- 

ful to Marianna Tsaplina for her assistance in preparing this 
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